Baltic Astronomy, vol. U, 000-000, 2005. 



LIMITATIONS OF THE HAMILTONIAN TREATMENT 
FOR COLLISIONLESS ASTROPHYSICAL ACCRETION 
FLOWS 



Vladimir I. Pariev 1,2 and Eric G. Blackman 3 

1 Physics Department, University of Wisconsin- Madison, 
1150 University Avenue, Madison, WI 53706, U.S.A. 

2 P.N. Lebedev Physical Institute, Leninsky Prospect 53, 
Moscow 119991, Russia 

3 Department of Physics and Astronomy, University of Rochester, 
Rochester, NY 14627, U.S.A. 

Received 2005 

Abstract. In the Hamiltonian treatment of purely mechanical sys- 
tems, the canonical and actual momentum of a particle are the same. 
In contrast, for a plasma of charged particles and electromagnetic fields, 
those two momenta are different. We show how this distinction is fun- 
damental in identifying the limitations of a recent attempt by Binney 
(2003) to rule out two-temperature collisionless astrophysical accretion 
flows from Hamiltonian theory. This illustrates the Hamiltonian method 
for astrophysical plasmas, its relation to the equations of motion, and 
its role in practical calculations. We also discuss how the complete 
Hamiltonian treatment of a plasma should couple the particle motion 
to a fully dynamical treatment of the electromagnetic fields. Our results 
stand independent from the discussion of Quataert (2003) who argued 
that time scale calculated in Binney (2003) is not the equipartition time 
as claimed. 
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1. INTRODUCTION 

An accretion disc of hot plasma orbiting a massive black hole 
and slowly spiralling in from the action of viscosity is widely consid- 
ered to be the main energy source for very luminous extragalactic 
objects such as quasars and other anomalously bright active galactic 
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nuclei (AGN) (Shapiro & Teukolsky 1983; Krolik 1999). Recently, 
the determination of the density and temperature of the gas in the 
accretion disc surrounding the central black hole of our own Milky 
Way Galaxy as well as discs around other nearby galactic centres has 
become possible from X-ray data taken by orbiting X-ray telescopes 
such as XMM Newton and the Chandra X-ray Observatory. From 
these measurements, the accretion rate, M, can be determined. The 
standard theory of geometrically thin accretion discs around black 
holes predicts that about 0.1 of Mc 2 is converted into escaping ra- 
diation. This estimate is consistent with the radiative efficiency of 
the integrated luminosity of quasars with the observed space density 
of supermassive black holes. However, a puzzling discrepancy has 
emerged: the measured luminosities of the central sources in some 
nearby galaxies are smaller than this standard estimate by 3 to 5 
orders of magnitude. 

This discrepancy led to the development of a newer, but also pop- 
ular, type of geometrically thick accretion disc models called "advec- 
tion dominated accretion flows" (ADAFs) ((e.g., Ichimaru 1977; Rees 
et al. 1982 Narayan & Yi 1995; Narayan, Mahadevan & Quataert 
1998) in which the gravitational binding energy of the accreting ma- 
terial is retained as internal energy within the hot plasma and ulti- 
mately crosses the event horizon of a black hole as the plasma falls 
in without radiating significantly. In the vicinity of the black hole 
horizon, the gravitational binding energy of the plasma is a fraction 
of Mc 2 and, therefore, the associated internal energy per particle in 
an ADAF is of order m p c 2 ~ 1 GeV. Because electrons are more 
mobile, they are the primary radiating particles. Thus, in order for 
such a weakly radiating accretion flow to exist, three main assump- 
tions of the ADAF model must be satisfied: (1) the internal energy 
dissipated in the accretion process via viscosity must go almost en- 
tirely into ions. (2) the heat transfer from ions to electrons must 
be slow enough, so that only a tiny fraction of the dissipated ther- 
mal energy received by the ions is transferred to electrons during the 
time it takes the gas to lose its angular momentum and fall onto the 
black hole. (3) the effective viscosity must be very high in order that 
the gas looses its angular momentum quickly and can indeed accrete 
faster than the ion-electron thermal coupling times. 

If the assumptions were true, an accretion flow with a given M 
could have a low enough number density and high enough tempera- 
ture that Coulomb collisions are inefficient in establishing equiparti- 
tion of energy between ions and electrons during the accretion time. 
In the absence of any other plasma process that could speed up the 
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ion-electron energy transfer, the electrons could remain at tempera- 
tures which are orders of magnitude lower than the ions in the black 
hole engine environments. Since electrons produce practically all of 
the radiation in a given disc, the luminosity rises steeply with increas- 
ing electron temperature. Thus keeping electrons at a much lower 
temperature achieves the main goal of ADAFs: substantially low- 
ering the luminosity for a given accretion rate (e.g., Ichimaru 1977; 
Narayan & Yi 1995; Quataert & Gruzinov 2000; Narayan, Igumen- 
shchev & Abramowicz 2000; Narayan 2002) compared to standard 
thin discs. ADAFs are thick discs because the heat dissipated by the 
accretion is stored as thermal energy of the ions, which puffs up the 
disc. 

Despite the important implications of the above assumptions if 
they were true, these assumptions have not been proven or disproven. 
Doing so requires understanding the subtle plasma physics of the 
interactions between ions and electrons with magnetic and electric 
fields. The assumptions have therefore been the subject of much 
deserved attention (Begelman & Chiueh 1988; Bitsnovatyi-Kogan 
& Lovelace 1997; Quataert 1998; Gruzinov 1998; Blackman 1999; 
Quataert & Gruzinov 2000). One central issue is whether or not 
collective long range interactions could be important for momentum 
transfer, energy dissipation and thermal equilibration processes oc- 
curring not just between pairs of particles (like Coulomb collisions) 
but in the whole volume of plasma, shortening the electron ion equi- 
libration time and ruling out the ADAFs. 

Since the basic microphysical processes of a plasma must involve 
known electromagnetic interactions between particles and fields, it 
is tempting to address the above ADAF assumptions starting from 
a very basic treatment of electromagnetic theory. This was recently 
pursued by Binney (2003), who presented a general argument that 
assumption (2) above is invalid. He used a Hamiltonian formalism 
to calculate the change of energy and angular momentum of a par- 
ticle moving in a general time-dependent electromagnetic field and 
averaged the result obtained over many particles to estimate the cor- 
responding rates for the plasma. He concluded that the ion-electron 
equipartition time t equ j is smaller than the characteristic angular mo- 
mentum loss time, t res (or residence time according to Binney's (2003) 
terminology). If this result were true, and the approach correct, it 
would rule out two-temperature accretion flows and close the whole 
ADAF chapter in astrophysics. 

Scientists working in plasma astrophysics usually start with writ- 
ing down Lorentz forces acting on particles and Maxwell equations 
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to analyse interactions in astrophysical plasmas. The idea of Binney 
(2003) to invoke general Hamiltonian analysis is novel and origi- 
nal. However, in this paper we follow along with the calculation 
of Binney (2003) and discover two fatal problems therein: (1) The 
calculation of Binney assumed that the particle angular momentum 
and the canonical angular momentum are the same, which leads to 
an incorrect and non-gauge invariant angular momentum equation. 
(2) The electric and magnetic fields are not included dynamically, 
as they must be for a plasma. We show that when the Hamiltonian 
analysis is performed correctly, no new conclusions can be made out 
of it that cannot be made by writing down usual Lorentz forces and 
Maxwell equations. Unfortunately, this exclude Hamiltonian analy- 
sis as a tool to resolve the problem of the existence of ADAF. We 
hope that a broader consequence of our consideration will help to 
elucidate the role and limitations of the Hamiltonian formalism in 
the context of such plasma dynamics problems. 



2. ANNOTATED DISCUSSION OF THE HAMILTONIAN 
APPROACH OF BINNEY (2003) 

2.1. The Rate of Change of the Hamiltonian 

Binney (2003) uses a Hamiltonian formalism to calculate the 
change of energy and angular momentum of a particle under the 
action of electromagnetic and gravitational fields. For a given time- 

— * 

dependent electromagnetic field characterised by the potentials A 
and and gravitational potential $, the Hamiltonian for the non- 
relativistic motion of a particle of mass m and charge q in Cartesian 
coordinates is (e.g., Landau & Lifshitz 1988b) 

ff= (P ~ g ^ )2 +# + m$, (1) 
2m 

where $ is assumed to be time-independent and axially symmetric, 
and the canonical momentum of the particle, p, is related to its 
velocity v by 

p = mv + qA, (2) 

where we set the speed of light c = 1 to match the notation of Binney 
(2003). Binney (2003) then calculates the time derivative of H. From 
general Hamiltonian theory it is known that the total time derivative 
of any quantity F along the trajectory of a particle is given by 

dF OF 

a-ar + wn (3) 
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where the Poisson bracket for position and momenta canonical vari- 
ables Si, Pi is {H, F} = ff-fp- — f^rffr with repeated indices summed, 
and dF/dt is the time derivative when Sj,pj are held fixed (Lan- 
dau & Lifshitz 1988a). Since {H,H} = 0, the rate of change of the 
Hamiltonian is 

dE _ dE p- qA dA dip 
dt dt ^ m dt ^ dt 

_ dA di) 

= - qv -^ + q m- (4) 



2.2. Hamiltonian is not the Particle Energy 

Binney (2003) then proceeds to identify dH/dt from equation (J3J) 
as a rate of change of the energy of a particle. He considers only 
the first term in the right hand side of expression (J3J to obtain the 
lower limit on the rate of energy transfer, and the upper limit on the 
equipartition time 

EE 

tequi ~ ]dEjdT\ < \ q v.dA/dt\ (5) 

These expressions need to be summed over the particles in some 
small volume to obtain the corresponding rates for the plasma as a 
whole. Assuming that individual positive and negative charges have 
the same charge magnitude \q\, we introduce the current density 
3 — \ a \ Hiwipt ~ ^m)i w h ere the sum is over ions and electrons in 
some local volume of plasma where is the velocity of the Ith 
positive charge, and w~ is the velocity of the mth negative charge. 
The result of separately summing up E and dE/ dt over the charges 
then leads to: 

jjtot 

tcqui < \Jd 3 xj -dA/dt\ (6) 

We must pause at this juncture. Binney's (2003) presentation of 
(JHJ) as a measure of the time rate of change of the particle energy is 
flawed because E is not in general a measure of the particle energy. 
Specifically, note that the Hamiltonian of Eqn. explicitly depends 
on time via A = A(f,t) and ij) = ip{r,t) and is not gauge invariant: 
the addition of a time derivative of some scalar function to ip changes 
E. Physical quantities such as the actual energy of a particle must 
be gauge invariant. It is not E that represents the kinetic energy 
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of a particle, but rather it is the first term in H that represents the 
particle kinetic energy. This term is explicitly gauge invariant (note 
that a gauge transformation also changes p) due to its dependence 

on A. It is because of these points, that the identification of the 
Hamiltonian with the energy of a particle and dH/dt as the work 
done on a particle by the electric field cannot be correct. Rather, it is 
the kinetic energy of a particle, K e = mv 2 /2, that must be estimated 
to determine equilibration times-the time derivative of the first term 
in Eq. ©. 

The equipartition time is not correctly obtained by (JHJ) but rather 
by t C qui = ^m e v 2 /\dK e /dt\, where dK e /dt is simply the time rate of 
change of the electron energy subject to acceleration by interaction 
with the electromagnetic fields. It is clear without detailed deriva- 
tion that dK e /dt is simply given by qv ■ E, the work done per unit 
time by the Lorentz force. Unlike in (jSJ) the correct expression for 
dK e /dt will have V<0 terms, which are particularly dominant when 
close interactions of electrons and ions occur (Coulomb collisions). 
Therefore, the correct estimate of t equ i is rather different from (jSJ) but 
is nothing new compared to the estimate of the conventional work 
done by the electric field. 

Binney (2003) further argues that in obtaining (JSJ) the qdip/dt 
terms of (HJ can be ignored because they cancel when averaging over 
particles with opposite charges in a quasi-neutral plasma. Actually, 
this cancelation is not guaranteed because a particle with a given 
charge in a plasma always attracts particles of the opposite charge 
within its Debeye sphere. As a result, the change in the potential at 
the location of a particle correlates with its charge and the cancella- 
tion of qdip/dt between species is not obvious. 

We have identified problems with © as presented by Binney 
(2003), but let us for the moment, continue to follow and dissect the 
arguments that he presents in comparing the equilibration and infall 
(angular momentum loss) times. We now discuss the latter and the 
pitfalls therein. 

2.3. Canonical Angular Momentum is not the Particle 
Angular Momentum 

To compare with the energy equilibration time ©, Binney (2003) 
proceeds to derive a similar time scale for the loss of angular mo- 
mentum of the plasma. If the angular momentum loss time were the 
longer of the two, he would rule out ADAFs because the electrons 
and protons would equilibrate before the plasma falls onto the black 
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hole. The fundamental problem with his subsequent calculation is 
that it relies on incorrectly identifying the canonical and particle 
angular momenta, as we now show. 

Specifically, Binney (2003) identifies the angular momentum of 
a particle L z with the 0-component of the canonical momentum p^. 
He uses Hamiltonian theory (e.g. Landau & Lifshitz 1988a) to write 

dL dH _ dA 

He then notices that the right hand side of expression (JJJ) differs 
from the right hand side of expression by replacing dt by dip and 
changing sign. Therefore, the analogous arguments and calculations 
leading to the estimate © are repeated to obtain an estimate of the 
"residence time" or the maximum time before the plasmas loses its 
angular momentum and falls into the black hole: 



tres ~ 



2 



dL**/dt\ \Jd 3 xj-dA/dcf>\ 



Dividing both sides of (JHJ) by both sides of (jHJ), approximating iJ tot /L* ot 
by Q<f, (acceptable when the plasma is substantially supported by Ke- 
plerian rotation as in the case of ADAFs), and approximating the 

ratio of integrals by the ratio of the corresponding derivatives of A 
one obtains 

^ < n/ffi* 1 . (9) 

tres * \dA/dt\ 

The ratio of derivatives of A can be approximated by the inverse 
frequency of a pattern of A propagating in direction. The lowest 
important frequency of the pattern in A is the Keplerian frequency 
Vtfj,. Then, the final estimate follows: 

tcqui < 1. (10) 



t 



res 



Taken at face value, this relation would imply that electrons receive 
a significant fraction of the thermal energy of ions before they have 
time to fall through the black hole event horizon, and ADAFs are 
impossible. We now pinpoint the key problem with the approach 
that led to this conclusion. 
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3. HAMILTONIAN EQUATIONS OF PLASMAS FROM 
FIRST PRINCIPLES 

The key problem with the result and calculation that leads to 
(II above is the mis-identification of the angular momentum of a 
particle L z with the 0-component of the canonical momentum p^ 
above Eq. (0). To show that these are not the same and the impor- 
tant consequences, we now rigourously derive dL z /dt starting from 
first principles. 

The forms of the Hamiltonian and canonical momenta depend 
on the choice of the coordinate system used to describe the particle 
motion. In cylindrical coordinates, expressions (fTJ) and (J2J) need to be 
modified. The general procedure to derive Hamiltonian equations is 
to start with a Lagrangian (which can be obtained from the covariant 
action). The Lagrangian is invariant under the transformations of 
space coordinates Sj. In cylindrical coordinates (si, S2, S3) =(r,<f>,z) 
and the Lagrangian C in the non-relativistic limit is (e.g., Landau & 
Lifshitz 1988b): 

£ = y if 2 + r 2 <p 2 + z 2 )+q (A r f + rA^ + A z z) 

—qip — m$, (11) 

where f = dr/dt, = d0/dt, and z = dz/dt along the trajectory of 
a particle. The conjugate canonical momenta are Pi = dC/dsi. In 
cylindrical coordinates 

p r = mr + qA r , p^ = mr 2 cf) + qrA^, p z = mz + qA z . (12) 

The actual particle angular momentum relative to the z-axis is L z = 
mr 2 (j). Eq. (fT2"j) allows us to relate L z to 0-component of conjugate 
momentum as 

L z =p lj) - qrA^. (13) 

Therefore, in electromagnetism L z 7^ p^ contrary to the statement of 
Binney (2003). As a result, his equation (JZJ) for dL z /dt is incorrect. 

To derive the correct expression for dL 2 /dt, we use the fact that 
the Hamiltonian H = J2i SidC/dsi — C. Then using (fTT|) and substi- 
tuting for Si from Eqs. (jT2"J) gives 



2m 



(p r - qA r ) 2 + {p^/r - qA (t> ) 2 + 
(p z - qA z ) 2 ] +qi/j + m$. (14) 
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This expression is different from the expression (JIJ for H in Cartesian 
coordinates by the factor 1/r multiplied with p^. Let us take d/dt 
along the particle trajectory of all terms in relation (|13|). One has 



dt 



q r—— + rep—— — h z 



{H, P(f> } 
dA z \ 



dH 



dip 



(15) 



where the curly bracket is again the Poisson bracket and we have 
used Eq. ()12|) . This result exactly reproduces the right hand side of 
equation (J7J) but the left hand side should be dp^/dt, not dL z /dt. 
There are extra terms to dL z /dt according to relation (J13j) : 



dL z 



dt dt 

dA^ , 

qr 



dA^ dA 4 
— qrA^ — qr—^ qr— — r 



dt 



dr 



dA d 



qr- 



dz 



-z. 



(16) 



Substituting — - from Eq. (J15|) into Eq. (|16j1 and collecting together 



dt 

terms with r and z we obtain 

dL z 



dt 



rz 



= q 
'idA 

r dd) 



. ldA r 
rr — - — 

\r dd) 

dAA 



dz 



r dr J 
/1# 
^ I r dd) 



+ 

dAt 
dt 



(17) 



When all components of potentials are combined in Eq. (J17|) . we are 
left with only components of magnetic and electric fields: 



-B z , 



ldA r 


1 d 


r dd) 


r dr 


ldA z 


dAtj, 


r dd) 


dz 


ldi/j 


dA^ 


r dd> 


dt 



Therefore, Eq. (IT?]) becomes 
dL, 



dt 



qrE^, + qr(v z B r - v r B z 



18) 
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This is nothing but the torque produced by the component of the 
Lorentz force qE + q(v x B) acting on a particle. Indeed such a 
result could be anticipated from the very beginning, without using 
the Hamiltonian or Lagrangian formalism for the derivation. In fact, 
equation (J7|) after Binney (2003) would already raise initial concerns 
from the fact that its left hand side dL z /dt, is a gauge invariant 
quantity, while the right hand side is not gauge invariant. 

The conclusion of Binney (2003) that we reproduced in the previ- 
ous section culminating in Eq. ()10|). crucially relies on the incorrect 
presence of the actual rather than canonical momentum on the left 
side of Eq. We have shown that when this equation is cor- 
rected, no new results come from the Hamiltonian formalism that do 
not already come from simply writing down the Lorentz force acting 
on a particle. The latter leads to classical estimates of equilibration 
times used in standard two-temperature accretion flow calculations 
(e.g. Narayan et al. 1998), and by contrast, the relation expressed 
in Eq. fjlO)) is simply invalid. 

4. INCLUSION OF THE ELECTROMAGNETIC FIELDS 

AS DYNAMICAL VARIABLES AND IMPLICATIONS FOR CAL- 
CULATING LOSS TIMES 

The energy transfer from protons to electrons in the turbulent ac- 
cretion plasma is mediated by electromagnetic fields (waves) excited 
by plasma instabilities. Particles and waves can exchange energy, 
(see for example Begelman & Chiueh (1988) in this accretion flow 
context). In order to quantitatively assess the effectiveness of the 
energy transfer between particles and waves, one needs to couple 
the dynamical evolution of the electromagnetic fields to dynamical 
evolution of the particle motion. This means that in general, the 
electromagnetic fields E and B (or A and ip) need to be treated as 
dynamical variables, not merely as background fields. This contrasts 
Binney (2003), who assumes that the electromagnetic field is fixed: 

A and i[) are considered as externally imposed background fields on 
the dynamical system of particles, not as dynamical variables them- 
selves. 

One might be tempted to extend the intended Hamiltonian ap- 
proach of Binney (2003) to include dynamical electromagnetic fields. 
In such an approach, the electromagnetic fields can be decomposed 
into the sum of normal modes and each such mode could be treated 
as a dynamical degree of freedom. One could then write a Hamil- 
tonian containing both particles and electromagnetic fields and use 
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it to calculate the time derivatives of the total energy and angular 
momentum of the system. But the result of such calculation is pre- 
determined: the total energy and angular momentum of all particles 
and electromagnetic fields (including radiation) must be conserved 
in time. Hamiltonian equations of motion written down for such a 
Hamiltonian will simply be Maxwell equations for the electromag- 
netic fields combined with the equations of particles motion under 
the action of Lorentz forces. Such a procedure would not lead to 
new results but simply lead us back to the conventional methods 
of analysing the equations of plasma dynamics which start with the 
equations of motion in the first place. 

5. CONCLUSIONS 

We have elucidated the role of the Hamiltonian formalism in the 
study of plasmas and its relation to the equations of motion. The 
latter is the usual starting point for the practical study of plasma 
dynamics in applications to laboratory and astrophysical plasmas. 
We have illustrated how the Hamiltonian formalism does not provide 
any more information than that which is contained in the equations 
of motion when it comes to the practical calculation of dynamical 
time scales of a system. 

We have illustrated the importance of understanding these basic 
derivations by means of an application to a very current topic in as- 
trophysics, namely, two-temperature, low luminosity accretion flows 
commonly used as an explanation for the otherwise mysterious qui- 
escent accretion engines at the centres of galaxies. In particular, we 
have discussed Binney's (2003) attempt to use general Hamiltonian 
methods to obtain a constraint on the ratio of ion-electron equipar- 
tition time to the angular momentum loss time of particles in these 
flows. If the former time scale were shorter than the latter, these 
accretion flows would be ruled out. 

We have shown that the mathematically correct Hamiltonian for- 
malism does not provide any new information for estimating the ion- 
electron equipartition time beyond conventional non-Hamiltonian 
approaches. This is revealed when one corrects the Hamiltonian 
approach of Binney (2003) by not equating the canonical angular 
momentum to the actual particle angular momentum. The revised 
calculation shows that the expression for angular momentum change 
of a particle used by Binney (2003) is incorrect (evidenced also by 
the fact that it is not gauge-invariant) and thus the subsequent con- 
clusion that the equilibration time between electrons and ions in 
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accretion plasma is always shorter than the accretion time is un- 
supported. Instead, if performed correctly, Hamiltonian expressions 
for the rate of change of angular momentum of a particle in elec- 
tromagnetic fields lead to the usual torque provided by the Lorentz 
force. 

We have also pointed out that to incorporate particle- wave inter- 
actions occurring in turbulent plasmas, one must treat the electro- 
magnetic fields as dynamical variables in the Hamiltonian formalism. 
Including these excitations will simply lead to a conserved Hamilto- 
nian, once again providing no new information beyond conventional 
plasma physics approaches. 

Finally, we note that Quataert (2003) also argued that the con- 
clusion of Binney (2003) is incorrect. Quataert (2003) argued that 
the time scale on which the energy of a particle changes due to the 
work by the electric field is not the time scale on which the true 
heating or change in entropy occurs. He mentions two examples 
where this difference is evident: First is the motion of a particle in 
a slowly varying magnetic field, with characteristic variation time 
much longer than where Q = eB/mc is the cyclotron frequency. 
After some time the magnetic field returns to its initial value every- 
where. In this case, t equ i calculated by the method of Binney (2003) 
[expression (jHJ) above] would be the characteristic variation time of 
the magnetic field. At the same time, in the absence of collisions, the 
energies of particles remain the same because of the conservation of 
the adiabatic invariant. His second example is an undamped Alfven 
wave. In this wave the energy is transferred periodically between 
fields and particles, but there is no net heating. As an extension 
of the argument about different time scales for adiabatic and dissi- 
pative energy changes, Quataert (2003) mentions that particles are 
heated at discrete wave-particle resonances, not explicitly accounted 
for in Binney (2003). 

Although these two examples of non-dissipative energy changes 
of particles in time variable magnetic fields are clear and correct, 
the statistical nature of the turbulence (presumably existing in any 
accretion flow due to the non-linear development of the magneto- 
rotational instability (MRI, e.g. Balbus & Hawley 1998)) does not 
allow one to conclude that all particle-turbulence energy exchange 
processes will occur as in Quataert's two examples. Therefore, by 
themselves, these arguments of Quataert (2003) do not disprove the 
derivation of Binney (2003). In particular, Binney argued that the 
rate of heating may be estimated from equation (2) in his paper [re- 
peated as equation (J2J) above in this paper] as follows: "Thus this 
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equation describes the mechanism by which equipartition is estab- 
lished between ions and electrons; the net direction of the energy 
flow is mandated by the general principles of statistical physics, and 
the rate of flow may be estimated from equation (2)." This state- 
ment does not contradict the specific energy transfer examples of 
Quataert described in our previous paragraph above. The reason 
is that it is not clear how statistically important the examples of 
Quataert are for a realistic accretion flow. We have found different 
and more fundamental reasons that rigourously disprove the results 
of Binney (2003). 
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